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SUMMARY 

A simple method is developed to solve plane-stress problems with 
axial symmetry in the strain-hardening range based on the defoliation 
theory of plasticity employing the finite -strain concept. The 
equations defining the problems are first reduced to two simultaneous 
nonlinear differential equations involving two dependent variables: 

(a) the octahedral shear strain, and (b) a parameter indicating the 
ratio of principal stresses. By multiplying the load and dividing the 
radius by an arbitrary constant, it is possible to solve these problems 
without iteration for any value of the modified load. The constant is 
determined later by the boundary condition. 

The method is applied to the cases of a circular membrane under 
pressure, a rotating disk without and with a hole, and an infinite 
plate with a circular hole. Two materials, Inconel X and 16-25-6, the 
octahedral shear stress -strain relations of which do not follow the 
power law, are used. Distributions of octahedral shear strain, as well 
as of principal stresses and strains, are obtained. These results are' 
compared with the results of the same problems in the elastic range. 

The variation of load with maximum octahedral shear strain- of the mem- 
ber is also investigated. 

The following results are obtained: 

1. The ratios of the principal stresses remain essentially con- 
stant during loading and consequently the deformation theory is appli-. 
cable to this group of problems. 

2. In the plastic deformation, the distributions of the principal 
strains, and of the octahedral shear strain, are less uniform than in 
the elastic case, although the distributions of the principal stresses 
are more uniform. The stress concentration factor around the hole is 
reduced with plastic deformation, but a high strain concentration 
factor occurs . 

3. The deformation that can be accepted by the member before 
failure depends mainly on the maximum octahedral shear strain of the 
material . 
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4. The added load that the member can sustain between the onset of 
yielding and failure depends mainly upon the octahedral shear stress - 
strain relation of the material. 


INTRODUCTION 

In the design of turbine rotors, it is desirable to know the 
detailed stress and strain distributions in the strain-hardening range 
and the amount of increase in load that can be sustained between the 
onset of yielding and failure. It is also desirable to know the effects 
of a notch or a hole in a turbine rotor or other machine members that 
are stressed in the strain-hardening range . If a member is thin, it can 
be analyzed on the basis of plane stress. For problems of this type, 

Nadai obtained solutions for ideally plastic material in the cases of 
the rotating disk, the thin plate with a hole, and the flat ring radially 
stressed (references 1 and 2). For the case of materials having strain- 
hardening characteristics, a solution of plane-stress problems has been 
obtained by Gleyzal for the circular membrane under pressure (reference 3). 
The concept of infinitesimal strain was used and the solution was obtained 
by an iterative procedure with a good first approximate solution. The 
plastic laws were always satisfied by using a chart given in reference 3. 

In reference 4, a trial -and -error method is given for rotating disk with 
very small plastic strain, in which the elastic stresses and strains are 
used as the first approximate values. Experimental investigation for the 
high-speed rotating disk is made in reference 5; distributions of plastic 
strains (logarithmic strains) for different types of disk are measured. 
Reference 6 experimentally investigates the burst characteristics of rotat- 
ing disks ; stress at the center of disk is calculated by assuming that the 
material behaves elastically at the burst speed; the average tangential 
stress along the radius at burst speed is also calculated. 

A simple method of solving plane -plastic -stress problems with axial 
symmetry employing the finite strain concept in the strain-hardening range 
and based on the deformation theory of Hencky and Nadai (references 7 to 
9), which is derived under the condition that the directions and the 
ratios of the principal stresses remain constant during loading, was 
developed at the NACA Lewis laboratory and^is presented herein. The 
equations of equilibrium, strain, and plastic law are reduced to two 
simultaneous nonlinear differential equations involving three variables, 
one independent and two dependent, that can be integrated numerically to 
any desired accuracy. These variables are the proportionate radial 
distance, the octahedral shear strain, and a parameter a that indicates 
the ratio of principal stresses. The magnitude of variation in calculated 
values of the parameter a with change in load directly indicates whether 
the deformation theory is applicable to the problem. 

The method developed is applied to: (l) a circular membrane under 
pressure, in order to compare results obtained by this method with that 
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obtained by Gleyzal (reference 3); (2) rotating disks without and with a 
circular central hole, in order to investigate plastic deformation in 
such disks and the effects of the hole; and (3) an infinite plate with a 
circular hole or a flat ring radially stressed, in order to investigate 
the effects of the hole in the strain-hardening range. 

In the investigation of (2) and (3), two materials, Inconel X and 
16-25-6, with different strain-hardening characteristics were used in 
order to determine the effect of the octahedral shear stress -strain 
curve on plastic deformation. The octahedral shear stress of these two 
materials is not a power function of the octahedral shear strain, so that 
more general information can be obtained. Distributions of stresses and 
strains of the same problems in the elastic range are also calculated for 
purposes of comparison. 


SYMBOLS 

The following symbols are used in this report: 
a radius of hole 

b outside radius of membrane, rotating disk, or flat ring 

c outside radius of plate, very large compared with radius a 

h instantaneous thickness of membrane, rotating disk, or plate 
h^ initial thickness 

k arbitrary constant 

p pressure on membrane' 

r radial coordinate 

s arc length 

u radial displacement 

v axial displacement 

z axial coordinate * 

a parameter indicating ratio of principal stresses 

y octahedral shear strain 

£ logarithmic strain 

* 

9 angular coordinate 

p mass per unit volume 

O normal stress , normal force per unit instantaneous area 
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t octahedral shear stress 

(0 angular velocity 


Subscripts : 
b at radius b 

c at radius c 

o at center, for case without hole or at radius a for case with 

concentric circular hole 

1,2,3 principal directions in general 

r,0,z principal directions in cylindrical coordinate system. 


STRESS -STRAIN RELATIONS IN PLASTIC DEFORMATION 

The deformation theory of plasticity for ideally plastic materials 
was developed by Hencky from the theory of Saint Tenant -Levy -Mise s for 
the cases in which the directions and the ratios of principal stresses 
remain constant during loading (reference 1 ). Nadai extended the theory 
to include materials having strain-hardening characteristics (refer- 
ences 8 and 9). The conditions for the deformation theory have been 
emphasized by Nadai (reference 9, p. 209), Ilyushin (references 10 and 
11), Prager (reference 12), and Drucker (reference 13). Experiments con- 
ducted by Lessells and MacGregor (reference 14), Osgood (reference 15), 
and others on thin tubes subjected to combined loads with the directions 
and the ratios of the principal stresses constant throughout the body and 
remaining constant, during loading show that good results can be expected 
from the deformation theory. 

In 'more recent experiments on thin tubes by Fraenkel (reference 16) 
and Davis and Parker (reference 17 ) , it has been shown that even with 
considerable variation of the ratio of principal stresses during loading, 
the strains obtained from the experiments were in good agreement with 
the strains predicted by use of the deformation theory. Further experi- 
mental investigation is needed to determine the extent to which the vari- 
ation of ratios of principal stresses is permissible . In case the vari- 
ation is small (approximately 10 percent over the strain-hardening 
range), the deformation theory can, however, be expected to give good 
results . 
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In the present problems with axial symmetry, the directions of 
the axes of principal stress remain fixed during loading and it seems 
that the ratios of principal strains and of principal stresses may 
also remain approximately constant. The deformation theory previously 
discussed is therefore used. The stress -strain relations are then as 
follows : 

e l + e 2 + c 3 = 0 (1) 


a l~ a 2 _ q 2~ g 5 _ <*5^1 
e l- e 2 € 2- e 3 e 3" e l 

T = t(7) (3) 

where 

t = | ^-ag) 2 .+ (a 2 -o 3 ) 2 + (a 3 ^ 1 ) 2 ] l/2 (4a) 

7 = 3 jj C l' € 2 )2 + (€ 2" € 3 )2 + («>) 

From equations (l) to (4b), the following relations are obtained: 

C 1 - | ? [?L - \ 

<2 ’ I 7 [°2 - I «W] 

‘3 ‘ 5 T [f = ' f 

For plane -3tress problems c 3 =0, so that 


^ 2 / 2 2^ 1//2 
T “ T ( °1 "°l a 2 +a 2 } 

(5a) 

[2 , 2 . 2, 1/2 
7 - 2 ^- (e x +£ie 2 -t€ 2 ) 


(5b) 
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and 

e l = 3 T ^ a l" . (6a) 

€ 2 “ 3 f (°2 ~ 2°1^ ( 6t ) 

£ 3 = |'T [" k = ' (£ l +€ 2 } (6c ) 

When Oj and c 2 are expressed in terms of e 1 and e 2 , there is 
obtained 


Oi = 2 I (2e 1+ € 2 )| 

a 2 = 2 I (2c. 2 +e 1 ) 


(7) 


EQUATIONS OF EQUILIBRIUM AND STRAIN INVOLVING DISPLACEMENTS 

Equations of equilibrium and equations of strain are derived for 
three plane -stress problems with axial symmetry. It is convenient to 
use cylindrical coordinates for these derivations; the principle 
directions 1, 2, and 3 in the preceding equations become radial, 
circumferential, and axial directions, respectively. Because a large 
deformation in the strain-hardening range will be considered, the 
concept that the change of dimension of an element is infitesimal 
compared with the original dimension of the element is not accurate 
enough. Hence, the finite -strain concept, which considers the instan- 
taneous dimension of the element, is used. (The equations of infin- 
itesimal strains will be given by considering them as special cases of 
finite strains.) The stress is then equal to the force divided by the 
instantaneous area and the strains are defined by the following equation: 


5(9 




where l* is the instantaneous length of a small element having the 

lJ 

original length of (Z^Jq and j = 1, 2, and 3. During plastic 
deformation, the plastic strains at a certain state depend on the .path 
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by which that state is reached. For the paths along which the ratios 
of principal stresses remain constant during loading, however, the 
octahedral shear stress -strain relation, the value of the octahedral 
shear strain, and the value of the principal strains are defined 
by the initial and final states (reference 15 and reference 9, p. 209); 
5(€j) is then an exact differential and 

£ 3- log =T^fe ■ ' (8) 

It should be noted that the condition under which equation (8) was 
obtained is also one of the conditions under which the deformation theory 
is derived; as long as the deformation theory is applicable, equation (8) 
can also be used . 


Circular Membrane under Pressure 

The membrane considered is so thin that bending stress can be 
neglected (reference 18, p. 576). Figure 1 shows the membrane clamped 
at the rim and subjected to a pressure p, and a small element defined 
by A 9 and As taken at radius r+u in the deformed state . In the 

undeformed state, the same element would be at radius r and defined 
by A 0 and Ar. The instantaneous thickness of the element and the 

stresses acting on the element are also shown in the figure . The two 
principal stresses are G r and Oq, and cp is the angle between O r 
and the original radial direction. 

Equations of equilibrium . - When all the forces acting on the 
element in the direction of are summed up, the following equation 

of equilibrium is obtained ; 

O r (r+u) h A0 - (° r +AO r ) |jr+u + A(r+u)J A0 (h + Ah) cos A<P + 

2 °q As (h + i Ah) sin — cos <¥> - p As (r+u) A0 sin ^ 

= 0 


When A (r+u) approaches- zero as a limit, the differential equation of 
equilibrium may be obtained: 

d(o h) 

(r+u) . <»> 
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A cap of the membrane bounded by radius r+u and the forces 
acting on it are shown in figure 2 . Summing up the forces in the 
z -direction yields 


or 


pit (r+u) 2 


a r ~ 2rt h(r+u) 



( 10 ) 


Equation of strain. - .Inasmuch as the element at radius r, defined 
by A© and Ar in the undeformed state, is moved to radius r+u, 
defined by A© and As, by the application of pressure, p (fig. 1), by 
, use of equation (8) the strains are 


!og e 


ds 

dr 


€ 0 = log, 


r+u 
>e r 

h 


S * log e h. 


Then 


^ j 1 * L , 


dw 1 
d(r+u) 


1/2 


• 0 r+u 


(11a) 

(lib) 


e z 

e 



(11c) 


Rotating Disk 

Equation of equilibrium . - A disk of radius b and thickness h, 
rotating about its axis with .angular speed od , and an element taken at 
radius r+u, defined by A© and A (r+u), . is shown in figure 3 with 
all the external forces acting on it. Summing up all forces acting on 
the element in the radial direction yields 



NACA TN 2217 


9 


a r (r+u) h A6 ■ (o r +Ag r ) Jjr+u + A(r+u)J A0 (h+Ah) 
3J e ^A(r+u)J (h + I Ah) sin ^ - 

r , 1 pn |(r+Ar) 2 - r 2 1 A 9 

<•> I r+u + g A(r+u)J 


V 


= 0 


When A(r+u) approaches zero as a limit, the following equation 
equilibrium is obtained: 


r+u 


<l(a r h) ' 22 

iH - h - h i 


r+u dr 
r d(r+u) 


Equation of strains . - The strains are 


e r = lo Se 


€ 9 = lo Se 


d(r+u) 

dr 

r+u 


z ~ l0S e h ± 


therefore 


d(r+u) 

dr 


£0 


r+u 


> e z _ h 


of 

( 12 ) 


(13a) 

(13b) 


(13c) 
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Infinite Plate with Circular Hole or Flat Ring Radially Stressed 

An infinite plate uniformly stressed in its plane in all directions 
and having a circular hole is shown in figure 4. The whole system is 
equivalent to a very large circular plate of radius c with a small 
concentric circular hole radially subjected to. the same uniform stress 
O on the outer boundary. The solution obtained in such a plate within 
any radius b can also be considered as a solution of a flat ring with 
outer radius b and inner radius a, that is, uniformly loaded at the 
outer boundary with the radial stress a-^ obtained in the plate 
solution. 

The equations for this case can be obtained in a manner similar to 
the two previous cases, or by simply setting dw/dr and w equal to 
zero in the case of the membrane, or setting u) equal to zero in the 
case of the rotating disk. 


EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY IN 


TERMS OF PRINCIPAL STRESSES AND STRAINS 
Circular Membrane under Pressure 

By combining the equations previously derived the following set of 
independent equations, which define the problem, are obtained: 


.17 
r 3 t 

<° r - ? °e’> 

(6a) 

1 7 
9 ~ 3 f 

(°e-| a r> 

(6b) 

1 7 
■z = 3 T 

[- 1 (' vw] 

(6c) 

r -4( 

r , 2 2NJ-/2 

^ € r +€ r € e +€ e J 

(5b) 

T = T(r) 


(3) 


* 
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a r _ d(r+u) A _dv 

dr \ d(r+u) 


12 


1/2 


•0 r+u 


h 


(r+u3 


d(o r h) 

d(r+u) 


= h(O e -Q n 


(11a) 

(lib) 


• (He) 



dw 

d(r+u) 


1_ 

' 2hO r 
p(r+u) 


2 


- 1 


(10) 


These equations are 10 independent relations of the 10 unknowns O , 

* X*' 

° 9 > c r> € e> € Z> y> V b, u, and w. 


•If equation (lib) is differentiated with respect to r and combined 
with equation (lla). 


r 


leg 

dr 


e ( £ r -%) 



~icT\ i/ 2 


dw 

a (r+u) 


J 


1 


Substituting equation (10) in equation (14) to eliminate 
following equation of compatibility: 


(14) 


w yields 



p(r+u) 



1/2 


1 


(15) 
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Equations (9) and (15) can "be simplified by using equations (11) to 
eliminate u and h, which results in 


Pa-°r) e (£ - £e) {l - 


rpe 


( e 0- € z) 


1 2j 1/2 


and 


r ~ - e ( € **- € e) J i - 


r dT - e y 

2hjO r 


2h i a r 

1 1/2 


- 1 


(16) 


(17) 


The ten equations defining this problem are now reduced to seven 
independent equations, (6a), (6b), (6c), (5b), (3), (16), and (17), 
with the seven unknowns a r , Cg, e r , €g, e z , T, and y. 


The solution of the problem is simplified by further reducing 

lj / 1 /> \ n 7n\ XT** . 


21 + a £ 

Mf) r Mi) 


= (o 0 -a r ) e (£r ' £0) l - 


r 

r- “| 

,2 

<! . 

pk r p ( e 0- c z) 
h^k 


V 

ro 

J 


E = e ( V C0) ^1 


k d /r 


(f) 


f 

pk r e ( € 0- e z) 

r 

1 

nj'k 

V 

t 

I 1 

L 2<7 r J 


l 


J 


- 1 


>(18) 


J 


where k is any arbitrary unknown constant with the dimension of 
length. By using the two parameters r/k and pk/hj^ it is then 
possible to solve the problem in a simple, direct way without use of 
the iteration method. This fact will be further discussed in the 
section METHODS OF NUMERICAL INTEGRATION. 
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Eotating Disk 

The set of equations that define this problem are: 

£ r = JT ( °r -l°e) 




17 r i 

3 T 2 


( a r + °e) 


- *$ (*r* ♦ 


T = t(7) 


J-t _ d(r+u) 


e = EtH. 
' r 

€ z h 
e = r— 


d(a r h) 

r+u dj^TT - *(a 0 -V 


,2 2. r+u dr 
saf r h ± — 


These equations are nine independent relations of the nine unknowns 
°r> °e> € r> € e> *Z> y> J > . h, and u. If equation (13b) is 

differentiated with respect to r and combined with equation (13a), 
the following compatibility equation is obtained: 

d€ e- (VV ' ' 

r d^ = e . " 1 (- 
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As in the case of the membrane, u and h can be eliminated 
from the equilibrium equation (12) by using equations (13), which 
yield 

r dir* + °r r dF - = ^" a r) e ^ r ^ “ P^r 2 e^" €z ^ (20) 


The nine equations defining this. problem are now reduced to seven 
independent, equations (6a), (6b), (6c), (5b), (3), (19), and (20), with 
seven unknowns a r , Oq, e r , e 0 , € z , t, and 7. 


The solution of the problem is made simpler by further reducing • 
equations (19) and (20) to the following forms: 





By using the parameters r/k and U)k instead of r and 00, a simple 
direct solution is possible for any arbitrary value of oak with k 
to be determined by the boundary condition. 


Infinite Plate with Circular Hole or Plat Bing Eadially Stressed 
The equations of equilibrium and compatibility for this case are : 


£ + O r L = (O r -OJ e (€ r -€ 0 ) ^ 

K 


*(§) Ml) 


(22) 


r d€ 6 = e ( e r- € e) 


k d /r 


(I) 




When equations (22) are combined with equations (6a), (6b), (6c), (5b), 
and (3), there are seven equations with seven unknowns. 
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EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY IN TERMS OF 


OCTAHEDRAL SHEAR STRAIN AND PARAMETER INDICATING 
RATIO OF PRINCIPAL STRESSES 

In the preceding section, displacements are eliminated from the 
equations, which result in seven equations involving the seven unknown 
quantities a r , Oq, £ r , eg, e z , • T , ' and y. The quantity e z can 

he expressed in terms of. e r and €g (from equation (l)). Two of the 
four unknowns a r , Og, e r , and €g may he eliminated hy using 

equations (6a) and (6b) or (7). The quantity T is a known function 
of 7 that is experimentally determined. The problem is then reduced 
to one involving three unknowns. Obtaining the solution of the result- 
ing equations is not, however, a simple matter. 

It is proposed that this difficulty can be avoided by using the 
following transformation: 


or 


a 0 -KT r = sin a 

c e~°r cos a 


O r { rjz sin a - cos a)\ 

a g =l ^i T <V* 8in a + 003 a v 


(23) 


Then O r and O q satisfy equation (5a). The octahedral shear stress 
T, a function of y, in the preceding equations varies- with r/k 
and also with loading. Such a transformation has been used for the 
Ideally plastic material (t = constant) by Nadai In the section "Yield 
ing in Thin Plate with Circular Hole or Flat Rings Radially Stressed" 
(reference 1, p. 189) and for a rotating disk (reference 2). From 
equations (6a), (6b), sand (23), the principal strains can be also 
expressed in terms of y and a; ■ 


£ = — —= (sin a - y 3 cos 

r 2 \j2 

£g = — ~ ; (sin a + a/H" cos 
“ 2/^2 


a) 


( 24 ) 
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The equations of equilibrium and compatihility for the three 
problems considered here are then obtained in terms of 7 and a in 
the following form: 


A 


A E-S2- + b tSL.. o 


k a/ r 


(I) Ml) 


Dlil.I *-&_ = F 

Ml) Ml) 


(25) 


J 

where the coefficients A, B, C, D, E, and F are functions of a, 
7, and r/k. For the circular membrane under pressure, from 
equation (18), 




A = 


B = 


7 cos a 


/ C 09 

3 sin a - cos a) — 

, r~ \/ 7 dT _ 7 sin 1 

(yj3 sin a - cos a)^ T dy ^ ) y 


C = 2 (cos a) e 


— 

(- | 7 cos a) 


Jf (,/j sin a t cos a) r 2 2 

1 ' ~ 2 ;f ^ On) 

6T 2 (^3 sin a - cos a) ' _ 


1/2 


D = (<^3 sin a - cos a) 7 
E = r (^ cos a + sin a) 


F = 2/^2 <1 - 


(-^j§7 cos a) 


f <V? 


sin a + cos a) 7 


1 - 


6T 2 (^3 sin a - cos a) 


rf/lT 2 


IS 




1/2 






J 


( 26 ) 
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For the rotating disk, from equation (21), 


A = ( \[5 cos a + sin a) - ( Jd sin a - cos a) - - 

V2 




B=(V3sl n a-=csa)(l^-m)i 


- \ cos 

C = 2(cos a) e ^ ' 2 


( - a /|>7 cos a ) ’ it- 2 ~ sin a 

W2 /-^|p (cok) 2 i(|) eV2 


> (27) 


D = (\/3 sin a - cos a) 7 
e = - ( Vi cos a + sin a ) 

r (~\li 7 COB a ) 

F = 2 \j2 1 - e\ ' 




For the infinite plate vith a circular hole, from equation (22), 
A = ( Vi cos, a, + sin a) - (\fs sin a - cos a) ^ a 


, rr . N / 7 dT 7 sin a\l 

B - (^3 sin a - cos a) gy ~^gr-)f 

f-/I 

(V 2 


7 cos a 


C = 2^cos 

D= (V5 sin a - cos a ) 7 
E - - (\/i " 008 a + 8 ^ n a ) 

F = 2^[l.e(^l 7CMa ) 


> (28) 




With these transformations, the solution of the problems is reduced 
to simply a numerical integration of the two simultaneous differential 
equations (equations (25)) involving the two unknowns 7 and a. 
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Furthermore, the parameter 7 , being the octahedral shear strain, . 
directly indicates the stage of plastic deformation at any point under 
any load, (in plastic problems, according to the deformation theory, the 
individual stress and strain distribution cannot give as clear a picture 
of the stage of plastic deformation as can the octahedral shear strain.) 
Also, the parameter a indicates the ratio of the principal stresses or 
strains. At any point, if a remains constant during loading, the ratio 
of principal stresses at that point remains fixed. The value of a 
obtained at each point in the calculation during loading directly indicates 
whether or not the deformation theory is applicable to the problem. 

The value of a is known at the boundaries or the center. From 
equations (23) and (24), in the case of a circular membrane under 
pressure , 

when r/b =0, 

O r = Oq 

a = | = 1.5708 - 

Ci 

when r/b = 1, 

€ 0 = 0 

a = fL « = 2.0944 
3 

In the case of a rotating disk without a hole, 
when r/b = 0, 


C r = Oq 
a = - = 1.5708 


0=0 


a = 2 . = 0.5236 . • 

6 


when r/b = 1, 



NACA TN 2217 


19 


In the case of a rotating disk with a hole, 
when r/a = 1 and r/b = 1 , 

O r = 0 

a = — = 0.5236 
6 

For the infinite plate with a circular hole, 
when r/a = 1 , 

q,. = 0 

a = — = 0.5236 
6 

when r/a approaches c/a or a value large compared with 1 , 

O r = a 0 

a = - = 1.5708 
2 


EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY FOR ^INFINITESIMAL STRAIN 

IN TERMS OF a AND 7 

The final forms of the equilibrium and compatibility equations for 
the. case of small strains are given in this section. The concept of 
infinitesimal strain is defined as follows: The change of dimensions, are 

small compared with the original dimensions, but are large enough so 
that the elastic strain can be neglected. The equations presented can 
be obtained either by direct derivation as was previously done or by 
reducing from the equations for finite strains through expanding the 

e f(a, 7 ) terms in series and neglecting the small terms. For inf ini - 
tesmal strain, the coefficients (functions of a and 7) A, B, C, 

D, E, and F in the preceding equations are denoted with a super- 
script prime in similar forms, but the coefficient (functions of a and 
7 ) are simpler than those for large strain. 
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A* £ J£L + B . E J2L 

k 4 (i) k *(§) 


D* + E' 21-32L- p- 


k d/r 


kd (l) 


For. the circular membrane under pressure, 


A' =1/3 cos a + sin a 

1 

B ' = ( y 3 sin a - cos a) i — 


C ' = 2 cos a 

D* = ( y? sin a - cos a )y 
E’ = -{\fz cos a + sin a) 

... # I 


F 1 = 2^3 7 cos a + 


Pk r 
k 


6 7(^3 sin a- - cos a,) 


For the rotating disk. 


A' =1/5 cos a + sin a 
B * = ( yi sin a - - cos a ) “ ^ 


C ' = 2 cos 


a-yfp te ) 2 (£) 2 i 


D * = ( Vs sin a - cos a)7 
E* = -(n/ 3 cos a + sin a) 
(cos a)7 
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For the infinite plate 

A' 

B' 

C* 

- D’ 
E' 
F » 


with a circular hole, 

= Vs cos a + sin a 

= {fJz sin a - cos a) i. 

T d.7 

= 2 cos a 

= ( /\/3 sin a - cos a) 

t 

= -(y3 cos a + sin a) 

= 2 ^ 3 (cos a)7 








(32) 


METHODS OF NUMERICAL INTEGRATION 

Two methods are developed to solve the differential eciations (25). 

In the first method, the differential equations are numerically integrated 
along r/k, which is considered the independent variable. In the second 
method, a is considered the independent variable . Because many terms 
in the equations are trigonometric functions of a, the use of a as 
the independent variable considerably reduces the work of computation. 


Numerical integration with r/k as independent variable . - 
Equations (25) can be written in the following forms: " — 


'N 


r da _ CE-FB 
k j/r\ AE-DB 


l (f) 


r dy _ FA -CD 
k j/r\ EA-BD 


l (D 


(33) 


J 


For the case of small strain, the terms A’, B’, C ', D*, E ' , 
and F* are used in equation (33) instead of A, B, C, D, E, and 


F, respectively. If at any point a and 7 are known. 


and 


da 

d(r/k) 

■ » f. can be calculated by equations (33) . At the boundaries or the 
<l(r/k) 

center, a is known, but 7 is to be determined by the load. Only 
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one value (unknown) of. 7 corresponding to a particular load exists 
on' each "boundary. It is therefore difficult to start the numerical 
integrations on the "boundary with the correct value of 7 corresponding 
to a given load. Also, in plastic problems covering the strain- 
hardening range, the method of superposition is invalid. Usually, a 
method of iteration is used to solve the problem (for example, refer- 
ence 3) . In the method presented herein, an arbitrary but unknown 
constant k has been introduced in equations (18), (21), and (22). For 
the cases considered, the terms in the equations that involve load are 


always multiplied by r, so that 


/pr\ : 


can be written as 


(if 


in equations (18) and (26) and (cor) 2 as (cok) 2 ( ^ in equations (21) 
and (27). Vlc/ 


The numerical integration can then be started at the inner boundary 
(or at the center if there is no circular hole at the center) by using 
the known values of a 0 , a desired value of y Q , and arbitrary value 

of |— \ for the membrane or of (cok) 2 for the rotating disk. The 
Vi/ 

numerical integrations can then be carried out, obtaining values of a 
and 7 at different values of r/k, until a progressively reaches 
the value that satisfies the other boundary condition. Because the 
value of r is known at the boundaries, the value of k can be 
determined for the selected value of 7 Q . The number of points and the 


formulas used in the calculation depend on the accuracy required (ref- 
erences 19 and 20). It lias been found that if the formula for evaluat- 
ing definite integrals is applied after using the forward integration 
formula (references 19 and 20), high accuracy can easily be obtained. 


The procedure used herein to obtain solutions is the same for each 
problem. Calculations are started from the inner boundary (or from the 
center if there is no circular hole at the center) with the known value 
of a Q , the desired value of y Q , and the arbitrary loading term. The 
parameter a Q is equal to it/2 at r/b = 0 for the membrane and for 


the solid rotating disk and is equal to rt/6 at r/a = 1 for the 
infinite plate with a circular hole and for the rotating disk with a 

hole. The arbitrary loading terms are [ — | and (cok) 2 for the 

\ h i/ 

membrane and the rotating disk, respectively. Then 


da 

( 

_ j v 


and 


[d$u| 0 ’ 


corresponding to a n and 


7 q at the inner boundary or the 
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center, are obtained from equation (33). The following formulas for 
forward integration are used to determine the first approximate values 
of a and 7 at the next point (a-^* and 7 *)• 



By substituting and 



7 2 * into equation (33), approximate values 
are obtained and the second approximate values 


of a,p and 
formulas : 


71 (a^** and 7 -j**) can be computed from the following . 



(34b) 



The values of or** and y~, ** are substituted into equation (33) again 

r” “T » — — » 


in order to calculated the values of 


By use of 


the following formulas for evaluating definite integrals, the values of 
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This procedure is applied to the next point, and so forth, until the 
value of a reaches the required value of at the outside boundary 

(aijj = 2/3 it at r/b = 1 for the membrane, 0 ^ = it/6 at r/b = 1 for 
the rotating disk, and = n/2 at r/a = c/a for the thin plate 
with a circular hole) . Inasmuch as 



the loading terms are determined as follows : 


For the membrane, 


For the rotating disk. 



(33a) 


(33b)* 


For the infinite plate with a circular hole, 


*c oh c (™) c - ° h l« 


-e. 


= Oh. ± e 


Jo 

2^2 


(33c) 


•or for the flat ring radially stressed at the outside diameter b, 

- -Zi. .(sin a. - Jz cos oO 

S * ("A h b - <°A V 


where t Q and t^ are the tension per unit original circumferential 
length at r = c and r = b, respectively. 
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Numerical Integration with a as independent variable. - 
Equations (33) can be written in tbe following forms: 


dr _ FA -CD 
da CE-FB 


ii). 


(35) 


AE-DB r 
da CE-BF k 




By using equations (26) to (28) and expanding 
the following equations are obtained: 


into a series, 


For the circular membrane, from equations (26), 
CE-BF = 2GjJ - |^2 /\/3HJ g^a,r, 7 - f- L (a,r) + 2 ^3 

AF-CD = 2 - 2BJ 7-2 \[ZL f 2 (a,r) j^a,r,|, 


4.^)1 


> (36) 


AE-BD = -l 2 - J 2 g 


( 7 dt\ 

V a ^7d7j 

For the rotating disk, from equations (27), 

CE-BF = -2BX. - 2 rfsEJ g^x,,7,I ^ fl(a,7) + L — f 3 (a,r) 

2 - 2 VSHL jl - f x ( a, 7 )J + J ^ f 3 (S7^ 7 

Mi) 


AF-CD = < 8H £ 


> (37) 


AE-BD. = -L 2 - J 2 g 
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For the infinite plate with a circular hole, from equations (28), 

CE ‘ ' 


-BP = - 2 HL - 2 ^3HJ g^a, 7 ,I ^( 0 , 7 ) 

AF-CD = fa E 2 - 2V3EL ji - f x (a,7^\ 7 


AF-3D = -L 2 - J 2 g 


( a ’ r 4 


dr 


J 


where 

G = sin a 
H = cos a 

J = y/z sin a - cos a 
V-V3 cos a + sin a 

= f &\ 2 


K l = 




=/| »<<*>’ 


% 

and. 

fl(a<,7) = 


^|(cos a ) 7 


1-e 


-J|(cos a) 7 




(38) 


= 1 - I ^|( cos a ) 7 + |(cos 2 a)7 2 - . . . 

-J|(cos a )7 
a, 7) * e 

= 1 - ,J|(cos a )7 + | X |(cos 2 a)7 2 - |^||( coa3 a )? 3 + • • • 
= 1 -#cos a)7 f x (a,7) 
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The symbols G, H, J, and L are trignometric functions of a 
only; K-j_ and Kg are constants during calculation. The symbols 

f±, f 2 > t and g are functions of a and 7; j is a function 

of a, y, and 

This method is used herein in the solution of an infinite plate 
with a circular hole. The procedure of numerical integration is 
similar to that used in the first method. The first four terms of the 

series of e ' f,/ are used; the accuracy of the result is the same 
as that in the first method, with a reduction of one half in computation. 

Both methods presented .herein are used to obtain the solutions for 
the given values of y Q . The purpose of. the present paper is to obtain 

solutions for the entire strain -hardening range and the methods developed 
are very convenient for this purpose. If, however, a solution for Only 
a particular value of loading is required, it can be obtained by inter- 
polating between values obtained from two or three solutions correspond- 
ing to loading near the specified value. 
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NUMERICAL EXAMPLES 

Membrane. - In order to compare the results obtained by the method 
developed herein to those obtained by Gleyzal (reference 3), one numeri- 
cal solution for infinitesimal strain is calculated by using the t( 7) 
curve of the tensile test in figure 1 of reference 3. Inasmuch as 
reference 3 states that: "For simplicity, strain will be taken to mean 

conventional. strain (ds-ds 0 )/ds 0 where ds and ds Q are final and 

initial arc length, respectively.", equations (29) and (30) for infini- 
tesimal strain are used. The calculation is started at r/k = 0.005. 

Values of a Q =' 1.5708, y Q = 0.0299, and pk/hj^ = 55,920 are used. 

Rotating disk . Numerical solutions for finite strain (equations (25) 
n.nfl (27)) are calculated. The t( 7) curves of two materials, Inconel X 
and 16-25-6, are plotted in figure 5(a)'. These data were supplied by 
W. F. Brown, Jr., H. Schwartzbart, and M. H. Jones. The same t( 7) curves 
are plotted on logarithmic coordinates in figure 5(b). These materials, 
Inconel X and 16-25-6, of which T is not a power function of, 7, were 
chosen so that more general information can be obtained. It should be 
mentioned that the given octahedral shear stress -strain curves (fig. 5) 
of these two materials have not been corrected for the triaxiality and 
nonnnifozm stress distribution introduced by necking, arid consequently 
do not represent the exact stress -strain relation after necking of these 
two materials. The solutions obtained from the T(y) curves of the 
tensile test after necking can, however, represent the solutions corre- 
sponding to materials having the exact i{y) curves shown in figure 5 and 
for simplicity the materials are herein still referred to as Inconel X 
and 16-25-6. 

. In each case, the calculation is started at r/k = 0.005, as in the 
case of a membrane. 

Three solutions are also obtained for a rotating disk with a central 
hole, using Inconel X. Calculations are started at r/a = 1. 
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All numerical examples are given in the following table : 


Solid rotating disk 


Material 


.IK p(o)b) 2 

1|3 . 

Inconel X 

0.04 

1 X 10 5 


.1152 

1 X 10 5 


.30 

1 X 10 5 

16-25-6 

0.04 

1 X 10 5 


.1152 

1 X 10 5 


.30 

2.5 X 10 5 


Rotating disk with central hole 


Material 


P(«a) 2 

Inconel X 

0.30 

1 X 10J 


.30 

2 x 10 5 


.30 

4 X 10 2 


Infinite plate with circular hole . - The calculations for this 
problem are carried out for the case in which C r = 0 at r/a = 1. 

The value of a Q at r/a = 1 is then 0.5236. (For other cases where 
a r is different from 0 at r/a = 1, the corresponding value of a 0 
should be used . ) The same materials as in the previous problem are 
considered. The numerical examples are: 


Material 

y 0 

Inconel X 

0.04 


.1152 


.1871 


.30 

16-25-6 

0.04 


.1871 


.30 
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RESULTS AND DISCUSSION 

The radial and circumferential stresses o r and Oq, respectively, 

obtained for the case of a circular membrane are plotted against r/b 
in figure 6. Two curves, taken from reference 3, corresponding to 
calculations for about the same pressure used in the present calculation, 
are included in the figure for comparison. In the present calculation, 
the t( 7) curve given in figure 1 of reference 3 and the same 
infinitesimal -strain definition based on the original dimension is used. 

In order to be consistent, the initial thickness h^ is also used in 
the calculation rather than the instantaneous .thickness h, which is 
used in reference 3. 

The variations of a with the radius for the rotating disk and 
for the infinite plate with a circular hole for different loads and 
materials are plotted in figures 7(a) and 7(b), respectively. The vari- 
ations of a with 7 0 (or loading) at various radii for the rotating 
disk and the infinite plate with a circular hole are plotted in 
figures 8(a) and 8(b), respectively. Similar curves for the ratio of the 
principal stresses O v /oq are shown in figures 9(a), 9(b), and 10. When 

figure 7 is compared with figures 9(a) and 9(b), it is seen that the 
variations of a with radius are very similar to the variations of 
O t /Oq with radius, although the relation between, a and O r /ag is not 

linear. 

Examples for a membrane with a large strain are not calculated 
herein, because the result of reference 3 is sufficient to give an approxi- 
mate variation of the ratios of principal stresses along the radius during 
loading, although the infinitesimal-strain concept is used. The variations 
of the ratio of principal stresses with radius for different loads, based 
on the values of O r and Oq given in figures 8 and 9 of reference 3, 

are calculated and plotted on figure 9(c). 

The values of C r are plotted against 0 Q at different radii 

under different loads for the rotating disk and the infinite *plate with 
a circular hole in figure 11. The heavy solid and dashed curves 
represent the values of a r and Oq at different radii for any given 
load and are called loading curves. The loading curve moves away from 
the origin with increasing load. The light solid and dotted lines con- 
necting the different loading curves at a given radius and extending to 
the origin represent the values of O r and Oq at different loads for 

any given radius and are called loading paths. Also shown in the fig- 
ures are the yielding surfaces, which are ellipses under the deformation 
theory. 
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A clear picture of the variation of the ratios of principal 
stresses in this group of problems with different loads and with 
different materials is given in figures 7 to 11. It is evident that 
the ratios of principal stresses remain essentially constant during 
loading. For this group of problems, the deformation theory is there- 
fore applicable and €j, defined by equation (8), represents the true 

strain and is determined by the final and initial states. 

The variations of y and y/y Q with radius are plotted in figures 12 
and 13, respectively, for the rotating disk and the infinite plate with 
circular hole. It is interesting to note that the curves in figure 13 for 
different loads for the same material are quite close. The curves for 
different materials on figures 7 and 9 are also close, but are not so 
dose in figure 13. . 

The- distributions of principal stresses and principal strains along 
the radius for the rotating disk and for the infinite plate with a 
circular hole are plotted in figures 14 and 15, respectively. For 
comparison, the variations of o 0 /(o 0 )b, € 0 /£ e ) t , 811(1 l/^b wlth 

radius for both the elastic and the plastic range are plotted in fig- 
ures 16 and 17 . (The equations for the elastic range are given in the 
appendix.) If only the stress distributions for the elastic and plastic 
cases are compared, it is seen that the stresses are more uniform in the 
plastic state; but if the distributions of the principal strains and the 
octahedral shear strain for the elastic and the plastic cases are 
compared, it is evident that a less -uniform strain distribution is 
obtained in the plastic state. It is of special interest in the case 
of the finite plate with a hole to note that with plastic formation, 
the stress (tangential stress) concentration factor around the hole is 
reduced, but instead there is a high concentration in principal strain 
and in octahedral shear strain. A similar conclusion regarding con- 
centration factor around a circular hole in a tension panel is obtained 
in references 21 and 22. 


The quantities O 0 /(a 0 ) o , a r /(o r ) 0 , e r /(€ r ) 0 , and e 0 /(e 0 ) o 

along the radius for a rotating disk and Oq/{Oq) q and e 0 / (e 0 ) o for 

an infinite plate with a circular hole are plotted in figures 18 and 19, 
respectively. The curves representing 0 T /{° r ) o > € r/^ e r^o # and - 

€ e/( € 0) o for different materials and different values of y Q are close 
together; but the curves of Oq/ (Oq) q are quite different for different 
materials, as well as for different values of 7 0 . 
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The relation "between load and maximum octahedral shear strain 7 
of the rotating disk and the infinite plate with a hole are plotted. in 
figure 20. The terms p(cob ) 2 and 1^/^ are designated the load for 

the rotating disk and the infinite plate with a hole, respectively. 

It is shown in these figures that the load increases considerably for 
Inconel X when the value of 7 ' increases from 0.04 to 0.30, whereas the 
load for 16-25-6 increases only slightly. 

In figures 7, 13, and 16 to 19, it is shown that in the case of the 
plate with a hole, the variations of a, ?/7 Q , € r /(€ r ) 0 > and Zq/(cq) q 

with radius are essentially independent of the value of 7 Q of the plate 
and the T (7) curve of the material, at least within the range of t, 7 
variation enclosed by those values of the two materials used. From these 
results, it can be seen that the deformation that can be accepted by the 
plate before failure depends mainly on the maximum octahedral shear strain 
(or ductility) of the material, which would not be true if the strain 
distributions were a function of the T-7 curve. In the case of the 
rotating disk, however, a slight effect of 7 Q and the t( 7 ) curve is 

apparent on the strains; this effect seems to be caused by the body- 
force term of the disk. 

The stress distribution that will determine the load the member can 
sustain is now considered. From figures 16 to 19, it can be seen that 
the variations of Oq/ (Oq) q with radius depend upon the t( 7 ) curve of 

the material and on the value of 7 0 of the member. From figure 20, it 
is also seen that the load depends on the t( 7 ) curve. It therefore 
follows that the added load that the member can sustain between the onset 
of yielding and failure depends on the t( 7 ) curve of the material. The 
octahedral shear (or effective) stress and strain curve of the material 
should therefore be used as a criterian in selecting a material for a 
certain member under a certain loading condition, because consideration 
of the maximum octahedral shear strain of the material alone (or ductility 
alone) is insufficient. 

The variations of a, 7 , a , Oq, e r , e@, and 7 / 7 ^ with 

radius of three rotating disks with a hole are shown in figure 21. These 
three disks have the values of ratios of outer and inner radius b/a 
equal to 5.32, 12.45, and 28.12, respectively. The tangential stress 
Oq, the tangential strain £ 0 , and the octahedral shear strain 7 are 

much less uniform then in the case of a solid rotating disk. The ratio 
of maximum and minimum octahedral shear strain 7 b is 'equal to 7.41 

for a disk with b/a = 5.32, equal to 11.75 for a disk with b/a = 12.45, 
and equal to 14.1 for a disk with b/a = 28.12; for a solid disk of the 
same material, the ratio 7 Q /y^ is about 5.3. 
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The values of the load, defined by p(wb)^, for rotating disks 
having different ratios of inner to outer radius a/h are represented 
by the 3olid curve in figure 22 . These disks were made of Inconel X 
and reach the same maxi mu m octahedral shear strain y Q = 0.3 at the 

inner radius of the disk. The dashed curve in the same figure is 
obtained by extending this solid curve toward a/b = 1, where the value 
of p((ob) 2 can be determined by considering a rotating ring with 
a/b — *1. The figure indicates approximately how the load p((x)b) 2 
varies with different disk's having different ratios of inner and outer 
radius and reaching the same maximum octahedral shear strain at the 
inner radius of the disk. The value of p(U>b) 2 for a.solid rotating 
disk made of Inconel X with y Q = 0.3 at the center of the disk is also 

indicated in the same figure . 

The preceding results and discussion were obtained for the plane - 
stress problems with axial symmetry in the strain-hardening range in 
which the elastic strains are negligible compared with the plastic 
strains. Whether these results and discussions are true for general two- 
dimensional or three-dimensional problems, or for the problems involving 
the region in which elastic strain is not small compared with plastic 
strain, or for partly plastic problems, can be determined only by a 
detailed. analysis of each case. 


CONCLUSIONS 

The results obtained in the cases of a membrane, a rotating disk 
without and with a hole, and an infinite plate with a hole in the strain- 
hardening range of two materials. Inconel X and 16-25-6, whose stress - 
strain relations do not follow the power law, show that: 

(1) The method developed not only accurately solved the, plane- 
plastic -stress problems with axial symmetry in a simple manner, but also 
gave a clear picture of the octahedral shear strain and the ratio of 
principal stresses during loading. 

(2) The ratio of the principal stresses in such cases remained 
essentially constant during loading and, consequently, the deformation 
theory is applicable to this group of problems. 

(3) The distributions .of principal strains, and octahedral shear 
strains, on the plastic state are less uniform than those' in elastic 

state, although the distributions of tangential stresses appear more uniform 
in the plastic state. The stress concentration factor around a hole is 
reduced in the plastic state, but instead there is a high concentration 
of principal strain and of the octahedral shear strain. 
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(4) The ratios of the strains along the radius of their maximum 
value are essentially independent of the value of the maximum octahedral 
shear strain of the plate and the octahedral shear stress -strain curve 
of the material. Hence, the deformation that can he accepted by the - 
plate before failure depends mainly on the maximum octahedral shear 
strain (or ductility) of the material. 

(5) The stress distributions depend on the octahedral shear stress- 
strain curve of the material. Hence, the added load that the member can 
sustain between the onset of yielding and failure depends mainly upon 
the octahedral shear (or effective) stress -strain curve in the strain- 
hardening range of the material. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronuatics, 
Cleveland, Ohio, February 28, 1950. 
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APPENDIX 


EQUATIONS FOE EOTATING DISK AND INFINITE PLATE WITH 
CIBCULAE HOLE IN ELASTIC EANGE 
Eotating Disk 

For a solid disk with the radial stress at the periphery (r = h) 
equal to zero, the principal stresses can he expressed in the following 
equations (reference 23, p. 68): 


= 

' r 8 


o„ = pa£(h 2 -r 2 ) 




o e = pofb 2 - pco2 r 2 


(39) 


where v is Poison's ratio. 

At r = h, 

(°9>h = \ P^ 2 (!- v ) 
Dividing equation (39) hy (cr^)^ yields 

°r _ 3+v 


"\ 


T 2(T3y 

°0 . 3+v 

T^b = 2X1^7 


1 - (f) 


1 - 


l+3v /r\' 
3+v 


(!) 


-‘J 


(39a) 


The stress -strain relations of plane -stress problems in the elastic 
range are: 


e r “ E ^r~ V£ ^ 


£ e = ^ ^e- u o r ) 


( 40 ) 
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where E is the modulus of elasticity in tension and compression. 


or 


Substituting equations (39) into equations (40) yields: 


e r = i (l-v)(3+v)(pa?b 2 ) 
^ (l-v)(,3+v)(pa£b 2 ) 


1 


1 


5(l+v) 

3+V 






«r 


3+V 

~2~ 


1- 


5(1+ v) /r\ 2 
3 +v \b/ 


^ e 0^h 


3+v 

2 


i-l ±1(L\ 


3+v 


(?) 


J 


(40a) 


(40b) 


The equations for the octahedral shear stress and strain given by 
equations (4a), (4b), and (5a) 1 can be applied to both the elastic and 
the plastic range, but equation (5b) is true only in the plastic range. 
The octahedral shear strain in the elastic range can be calculated by 
equation (4b) or simply by using the following equation: 


= 5%a T= ih±a|(o r 3.a r0e+ a 8 b l/2 


( 41 ) 


Substitute equations (39) in equations (41) to obtain: 


7 = 


or 


12E 


(pwV)- (3+V) 2 - 4(1+ V) (3+ v) + (7+2v+7v) (?) 


1/2 


7b = 2T 


— j - |”( 3 + V ) 2 - 4 ( 1 + 1 *( 3 + V ) (|) 2 + ( 7 + 2 U + 7 V 2 ) (|) 4 


1/2 


(41a) 


(41b) 
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and 

v = 0.29 for Inconel X (reference 24) 
v = 0.286 for 16-25-6 (reference 25) 


Infinite Plate With Circular Hole 

For a uniformly loaded infinite plate with a circular hole, the 
principal stresses are (reference 23, p. 56): 


A 

■N 



a r = -pr 

+ 2C 

V 

► 

(42) 

■ °e = “T 

+ 20 



r 

J 



where A and C are arbitrary constants. For the 

case considered 

herein, the boundary conditions are: 




, °r = 0 

at r = 

a 


°r = (a A 

at r = 

b 

• 

These boundary conditions are used to determine 

the arbitrary , 

constants A and C, which yield 




2 > 



<°r>b 

( I )- 1 



0 r = „2 




H 

1 

(?) 

/rx 2 

> 

(42a) 

' /-P _ - L '_ . 

(?) +1 



°0 .2 

, „ .2 



HD 

(?) , 
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Substituting equations (42a) into equation (41) yields 
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Figure 5. - True octahedral shear stress-strain curves. 
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Figure 7. - Variations of parameter a with proportionate radius for Inconel X and 16-25-6. 
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Figure 7. - Concluded. Variations of parameter a with proportionate radius for Inconel X and 16-25-6. 
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Proportionate radial distance, r/b 
(a) Rotating disk. 

Figure 9. - Variations of ratio of principal stresses with proportionate radial distance. 
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Proportionate radial distance, r/a 
(b) Infinite plate with hole. 

Figure 9. - Continued. Variations of ratio of principal stresses with proportionate radial distance. 
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Figure ||. - Loading curves, loading paths, and yielding surfaces. 
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(t>) Infinite plate 
Figure II. - Concluded. Loading curves, 


with circular hole. 

loading paths, and yielding surfaces. 
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(a) Rotating disk. 


Figure 12. - Variation of octahedral shear strain 
with proportionate radial distance. 
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(a) Rotating disk. 


Figure 13. - Variation of ratio of octahedral shear strain 
with proportionate radial distance. 



Proportionate radial distance, r/a 
(b) Infinite plate with hole. 

Figure 13. - Concluded. Variation of ratio of octahedral shear 
strain with proportionate radial distance. 
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(a) Rotating disk. 


Figure 15. - Variation of radial and circumferential strains 
€ r and e@ with proportionate radial distance. 
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(b) Infinite plate with holes. 


Figure 15. - Concluded. Variation of radial and circum- 
ferential strain e r and c g with proportionate 
radial distance. 
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(a) Variation of ratio of tangential stress to minimum 
tangential stress at outer radius vith proportionate 
radial distance. 



Proportionate radial distance, r/b 

(b) Variation of ratio of tangential strain to minimum 
tangential strain at outer radius vith proportionate 
radial distance. 

I 

Figure 16. - Comparison of results obtained in elastic 
and plastic range for rotating disk. 
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(c) Variation of ratio of octahedral shear strain to minimum 
octahedral shear strain at outer radius with proportionate 
radial distance. 


Figure 16. - Concluded. Comparison of results obtained in 
elastic and plastic range for rotating disk. 
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(a) Variation of ratio of principal stresses to tangential 
stress at outer radius with proportionate radial 
distance. 


Figure 17. - Comparison of results obtained in elastic and 
plastic range for infinite plate with circular hole. 
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(b) Variation of ratio of principal strains to tangential 
strain at outer radius with proportionate radial 
distance. 


Figure 17, - Continued, Comparison of results obtained in 
elastic and plastic range for infinite plate with circular 
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(c) ( Variation of ratio of octahedral shear strain to 
octahedral shear strain at outer radius with propor- 
tionate radial distance. 


Figure 17. - Concluded. Comparison of results obtained in 
elastic and plastic range for infinite plate with circular 
hole. 
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(a) Variation of ratio of principal stresses to stress at 
center of disk with proportionate radial distance. 


Figure 18. - Comparison of stress and strain distributions 
for different materials and for different maximum octahe- 
dral shear strain of rotating disk. 
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Proportionate radial' distance, r/b 

(b) Variation of ratios of principal strains to strain at 
center of disk with proportionate radial distance. 

Figure 18. - Concluded. Comparison of stress and strain 
distributions for different materials and for different 
maximum octahedral shear strain of rotating disk. 
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Proportionate radial distance, r/a 

(a) Variation of ratio of tangential stress to tangential 
stress at inner radius with proportionate radial distance. 

Figure 19. - Comparison of stress and strain distributions 
for different materials and for different maximum 
octahedral shear strain of infinite plate with circular 
ho 
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Maximum octahedral shear strain, y 0 


(h) Flat ring radially stressed. 

Figure 20. - Relation between load and mavlnmm 
octahedral shear strain. 
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Proportionate radial distance, r/b 
(b) Variation of octahedral shear strain with 


proportionate radial distance. 

Figure 21. - Continued. Rotating disk with hole, 
Inconel X, y 0 = 0.30. 
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(d) Variation of principal strains with 
proportionate radial distance. 


Figure 21. - Continued. Rotating disk with hole 
Inconel X, y Q = 0.30. 
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Ratio of inner and outer radii, a/b 


Figure 22. - Variation of load (function of speed) with ratio 
of inside and outside radii of rotating disk with hole. 
Inconel X; y 0 = 0.3. 
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